Probing Localization in Absorbing Systems via Loschmidt Echos 
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We measure Anderson localization in quasi-one-dimensional waveguides in the presence of ab- 
sorption by analyzing the echo dynamics due to small perturbations. We specifically show that the 
inverse participation number of localized modes dictates the decay of the Loschmidt echo, differing 
from the Gaussian decay expected for diffusive or chaotic systems. Our theory, based on a random 
matrix modeling, agrees perfectly with scattering echo measurements on a quasi one-dimensional 
microwave cavity filled with randomly distributed scatterers. 
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The propagation of waves through complex media is 
an interdisciplinary problem that addresses areas as di- 
verse as light propagation in fog or clouds, to electronic 
and atomic-matter waves used to transmit energy and 
information. Despite this diversity, the wave character 
of these systems provides a common framework for un- 
derstanding their transport properties and often leads to 
new applications. One such characteristic is a wave inter- 
ference phenomenon. Its existence results in a complete 
halt of wave propagation in random media which can be 
achieved by increasing the randomness of the medium. 
This phenomenon was predicted fifty years ago in the 
framework of quantum (electronic) waves by Anderson 
[l[ and since then has developed as a field of its own. 

Despite the enormous research efforts by various 
groups in measuring Anderson localization, it took nearly 
40 years to observe localization phenomena beyond any 
doubt. A decisive step towards this direction was done by 
optics and microwave experiments which allow a detailed 
study of the Anderson localization, undisturbed by inter- 
actions or other effects which characterize electron prop- 
agations. First experiments showing photon localization 
[2| had the problem of separating localization from ab- 
sorption, which can be another source of exponential de- 
cay of a propagating electromagnetic wave. A solution 
to this problem was given by Chabanov et. al in Ref. [3|, 
where they proposed to study the relative size of fluc- 
tuations of certain transmission quantities. They found 
clear evidence of localization in a quasi-one-dimensional 
(ID) microwave waveguide with randomly distributed di- 
electric or metallic spheres |3| (see also Fig.rjJ). 

This approach, although quite powerful, does not allow 
the view of transport from a dynamical perspective, nor 
makes a direct contact with the original ideas of Ander- 
son theory, which suggests probing localization by means 
of the sensitivity of the system properties against small 
perturbations [J]. This approach led us in recent years 
to focus on new measures that efficiently probe the com- 
plexity of quantum time evolution. One such measure 
is the so-called Loschmidt echo (LE), or fidelity, which 
probes the sensitivity of quantum dynamics to external 
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FIG. 1: Left Inset: Experimental apparatus - the top plate 
with two mounted antennae has been removed to show the 
waveguide with scatterers. Main Figure: The variance of the 
normalized transmission intensity f = S21I 2 / (|<S2i| 2 ). Ver- 
tical dashed lines are the mode cut-off frequencies where the 
number between is the number of open modes. The horizontal 
dashed line is the critical localization threshold of 7/3. The 
shaded region is a localized frequency window of 6.0—7.5 GHz. 
Right Inset: Normalized transmission distribution V(T) for 
the localized frequency window of 6.0 — 7.5 GHz. The red 
dashed line is a fit of the core to a log-normal distribution. 



perturbations (for review, see Ref. Q). The recent lit- 
erature on the subject is vast and ranges from atomic 
physics @ , microwaves Q , elastic waves to quan tum 



information [9(, and quantum chaos [lOl 
Formally, the LE F x (t), is defined as [3] 
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where f x (t) is the fidelity amplitude, H\ — Hq + XV is a 
one-parameter family of Hamiltonians, Hq is the unper- 
turbed Hamiltonian, AV^, where (|V^ TO | ) = 1, represents 
a perturbation of strength A, and |^o) is an initial state. 

Fidelity in its original definition, Eq. |T]) , is hardly ac- 
cessible to any experiment where the information about 
the system's state is based on the measurements of cer- 
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tain observables; the most popular being the scattering 
matrix itself. Therefore, the notion of scattering fidelity 
had been introduced as an alternative to Eq. (TP} [15| . In 
fact, it was shown that under certain conditions the scat- 
tering fidelity coincides with the standard fidelity [151 ]. 

We present here the first measurements of fidelity of lo- 
calized waves. Using a quasi-lD disordered cavity in the 
localized regime (see left inset of Fig. [I]), we investigate 
the fidelity decay of microwave radiation, due to small 
perturbations A in the form of boundary displacements 
of the sample, and find deviations from a Gaussian decay 
expected in a frequency interval associated with extended 
waves. Using a banded random matrix theory (RMT) 
modeling, we find that the fidelity amplitude decays as: 



f(t) ~ (ai) 2 csch 2 (crf); a = \y/l.5I 2 



(2) 



where I2 = J |"0(r)| 4 <ir is the inverse participation num- 
ber (IPN). Using a scaling analysis of a with respect to 
A we extract I2 and measure the localization properties 
of the sample, even if absorption is present. Our theoret- 
ical results, Eq. (|2|), are confirmed by our experimental 
measurements of scattering fidelity. 

The scattering fidelity amplitude is defined as 

fa b {t) = (s^(t)s° ab (t)) J ^(\s^ b (t)\ 2 )(\s° ab (t)f) (3) 

where S a b(t) is one component of the scattering matrix 
in the time domain, while the indices a and b refer to 
the antennae involved. In microwave studies the S ab are 
directly accessible from transmission or (for a = b) re- 
flection measurements in the frequency domain. The 
corresponding quantities in the time domain are then 
obtained by Fourier transforms. The super-indices A, 
indicate scattering matrix elements corresponding to the 
perturbed and unperturbed system, respectively. It is 
important to point out that the denominator in Eq. Q 
renormalizes the fidelity decay due to absorption, thus 
allowing us to trace out localization phenomena 15 1. 

The experiment has been performed in a quasi-lD rect- 
angular waveguide (height 8 mm, width 100 mm, length 
1190 mm) containing 186 randomly distributed brass 
cylinders of radius 5 mm (see left inset of Fig. [IJ . One 
end of the waveguide holds a fixed reflecting metallic 
wall, while at the other end there is a moveable reflecting 
metallic wall, which can be adjusted by means of a step 
motor. One antenna was placed close to the moving wall, 
while another one was placed in the center of the scatter- 
ing arrangement. Measurements have been performed in 
the frequency range 3 to 12 GHz. The cut-off frequency 
for the lowest mode is 1.5 GHz, while at the upper limit 
of the studied frequency range there are 7 propagating 
modes. The reflection amplitude at the center antenna 
(S22) and the transmission amplitudes between the two 
antennae (S21) have been measured for different wall po- 
sitions. The reflection amplitude S22 is used in Eq. ([3]), 
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FIG. 2: Experimental scattering fidelity in the frequency win- 
dow of 6 to 7.5 GHz (localized). The dots correspond to the 
experimental scattering fidelity, Eq. ©, for a wall shift of 
0.8 mm. The dashed line is a best fit to Eq. Q, while the 
solid line is a best fit to Eq. Both fits were done for the 
region t <tn- Eq. ((2)1 can be seen as the better fit. 



and the transmission amplitudes £21 are used in trans- 
missive studies below. The wall shift has been performed 
in steps of Sw = 0.2 mm up to a total shift of 18 mm. In 
addition, an ensemble average over 15 different realiza- 
tions of scatterer positions has been performed. 

In order to get confidence that our analysis is per- 
formed within the appropriate frequency window where 
localization is present, we first investigate the vari- 
ance o~~~ of the normalized transmission intensity T = 

\S2i | 2 / (\S2i | 2 ) [H- Since our experiment does not probe 
the total transmission but just one component of the 
scattering matrix, we expect localization whenever cr;~ 
exceeds the critical value of 7/3 Q. We find (see Fig.Qj) 
that this condition is satisfied approximately in the fre- 
quency window 5.5 — 9 GHz. Above 9 Ghz the waveguide 
modes are delocalized, while below 5.5 GHz the values of 
the variances are error prone, as £21 is below the pre- 
cision of the vector network analyzer (IS21I < 10~ 6 ). 
In the delocalized regime random matrix theory pre- 
dictions are applicable [lfij . yielding a value of a 1 ^ — 
(27V+l) 2 /[iV(2iV+3)]-l, where N is the number of open 
channels. In the limit N — > 00, the variance approaches 



the value cr~, = 1, in agreement with our experimental 
data for the high frequency regime (see Fig. [1} . 

We have also investigated the whole normalized trans- 
mission distribution V{T). As expected [l7[, we find 
a transition from a Rayleigh-likc behavior (applicable 
in the delocalized regime) to a broader distribution ap- 
proaching a log-normal behavior deep in the localized 
regime (see right inset of Fig. [1] showing V(T) in the lo- 
calized frequency window of 6 — 7.5 GHz.). 

Next, we measure the scattering fidelity in the localized 
frequency window (see previous analysis), and compare 
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the experimental data with the LE decay law found for 
chaotic/diffusive systems [5, 7]. The latter reads 



F(t) 



-\ 2 C{t). 



C(t) = t 2 +t- / dr / dT'b 2 (T r ), (4) 



which for small perturbations can be approximated as 
F(t) ~ cxp(— (Ai) 2 ). Above, 62 is the two-point form 
factor for the Gaussian orthogonal ensemble [5( . Using A 
as a fitting parameter, we have attempted to fit the ex- 
perimental data with Eq. (H|) . We found that the overall 
agreement is poor (see Fig. [5]). Further analysis (see be- 
low) confirms that Eq. (0| is inapplicable in the localized 
regime. On the contrary, when fitting the experimental 
data with Eq. ((2]) we get an excellent agreement. 

The analytical calculation of f(t) relies on the relation 
between scattering fidelity and LE in the weak coupling 
regime where our experiment operates. The LE was fur- 
ther evaluated using a RMT modeling for Hq and V. 
For diffusive or chaotic cavities, Hq and V are modeled 
by matrices drawn from a Gaussian Orthogonal Ensem- 
ble (GOE). Anderson localization in quasi- ID disordered 
systems, on the other hand, is modeled by sharp banded 
GOE matrices with bandwidth b for Hq and V 18, 1^ . 
In this model, the localization length scales as ~ b 2 . 
Therefore localization sets in for matrices of rank L 3> b 2 . 

Expanding the initial preparation as |-0o) = 



where Hq\ 



E°\4>°), it is found that the 



fidelity amplitude defined in Eq. |(Tj) can be written as 
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where H\\^) = E*\(f>n)- Averaging over disordered re- 
alizations one further gets that 



/(<) = £mn exp[-i(£ft - E°)t] 



(6) 



where a = J 1 - En<z IV'oMI 2 , ipo{n) are the components 
of the initial preparation \i/jq) m the position (Wannier) 
basis, and C mn = \ ((/>£, ]<ftS ) | 2 is the Local Density of 
States (LDoS) kernel [201]. In order to derive Eq. ([6|) 
we have assumed statistical independence between the 
eigenf unctions and eigenvalues of our Hamiltonian, while 
localization enforces the following contraction rule 



a* n a k = a 5 n ^. 



(7) 



For small enough perturbations, such that only nearby 
(on the order of mean level spacing) energy levels are 
mixed, one can approximate C mn in Eq. ^ with a delta- 
function. The fidelity amplitude is then written as 



dvV(v)exp(-iv\t) (8) 




FIG. 3: Numerical fidelity, Eq. Q, for Hamiltonians modeled 
by banded random matrices. The parameters are A = 10~ 3 , 
L = 5000, b — 10 (upper curve), and b = 3 (lower curve). 
Circles are the numerical results. Dashed lines are best fits to 
Eq. @. Solid lines are best fits to Eq. P). which again show a 
better fit. The inset shows the variance extracted from the fit 
of Eq. ([5} , plotted against the root of the inverse participation 
number. The observed linear relation confirms the validity of 
the theoretical calculations. 



where v n = {E*-E° n )/\ are the so-called level velocities, 
and V(v) is their corresponding distribution. The latter 
has been calculated in Ref. 2l| , and it was found that in 
the localized regime it is given by the expression 



P{t]) 



7T TTT] COth(7T?7/v / 6) — V6 
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sinh 2 (-7T 77/^6) 



(9) 



where r\ = v/a v . For localized level velocities, the vari- 
ance is a v — \fh, [22]. Substituting the above distri- 
bution in Eq. ([8]) we finally get the expression given in 
Eq. ([2|). The latter is checked numerically for two differ- 
ent bandwidths b = 10, 3 and for a system size L = 5000. 
Our numerical results are reported in Fig. [31 together 
with Eqs. ([2l QJ. To confirm further the validity of our 
calculations, we have fit the decay of LE for various band- 
widths b, with Eq. ©. From the fit we have extracted 
a v , which we have plotted against the IPN J 2 , in the in- 
set of Fig. [3] The observed linear behavior gives further 
confirmation to our theoretical calculations. 

We then analyze the decay of scattering fidelity due to 
small displacements of one wall of the disordered quasi- 
1D waveguide shown in the inset of Fig. [T] The shift of 
the wall can be mapped onto an effective Hamiltonian 
H\ with matrix elements [151 ] 



(Hx w ) nm = -wJ v±[Mv)W±[1>M]dv (10) 



where w and L are the shift and length of the moving 
wall, respectively, and Vi^n and V j_ip m are the normal 
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FIG. 4: Within the localized frequency window of 6 — 7.5 GHz, 
the experimental fitting parameters - a (squares) from Eq. ([2]) 
and A (circles) from Eq. Q - are plotted against the rescaled 
wall shift w/8w. Straight lines indicate best fit to power laws, 
a ~ w 0,92 , A ~ w . Inset: Same as main figure, but for the 
delocalized frequency window of 10.5 — 12 GHz. As opposed 
to the main figure, here A ~ w 1,a , in agreement with Eq. Q. 



derivatives of the wave functions at the wall. Thus a\ w cx 
w 2 . The proportionality constant is (2Lk 3 /Sir 3 ), for the 
case of chaotic cavities in the semiclassical limit [l^]. In 
any case, we finally get that X w ~ w. 

Next we proceed by fitting our experimental data on 
scattering fidelity with Eqs. ([2l HJ) where a, A are used as 
fitting parameters. We have extracted a, A, for various 
wall shifts w and plotted them versus a rescaled shift, 
w/5w. The results are summarized in Fig. [4j We find 
that in the frequency window 6 — 7.5 GHz (where An- 
derson localization is present), the best fit with Eq. (|4|) 
gives A ~ AJ,, with 7 = 1.9 ± 0.05. This result violates 
the theoretical expectation 7 w 1 and constitutes a direct 
confirmation that the RMT result Eq. (0| is not applica- 
ble in the localized regime. At the same time, the best fit 
of the experimental scattering fidelity with the prediction 
of the banded RMT modeling ||T]) gives that a ~ AJ, with 
7 ~ 0.92 ± 0.05, in agreement with our theory. The ex- 
tracted slope a/X w ~ \fli can be used as an estimation 
for the localization properties of our sample. Within the 
delocalized frequency window 10.5 — 12 GHz, a fit with 
Eq. (|4]) works perfectly well with 7 w 1.0 ± 0.05, in nice 
agreement with theory (see inset of Fig. 2]). 

In conclusion, using echo dynamics, we were able to iso- 
late absorption phenomena and identify unambiguously 
traces of localization in random media. Our theoretical 
calculations, based on a RMT modeling, indicated that 
the IPN of localized modes dictates the behavior of LE 
which follows a novel decay law; being experimentally 
distinguishable from the Gaussian decay observed in dif- 
fusive/chaotic systems. Our experimental measurements 
with disordered waveguides confirm the theoretical ex- 



pectations, thus suggesting fidelity as a reliable tool to 
investigate localization in the presence of absorption. 
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